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Abstract 

Another stabilizer term is used in the classical Hamiltonian of the 
Skyrme Model that permits in a much simple way the generalization of 
the higher-order terms in the pion derivative field. Improved numerical 
results are obtained. 
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It is well established that the Skyrme Model reproduces with relative suc- 
cess most of the static properties of Nucleon (within approximately ~ 30%). 
The idea consists of treating baryons as soliton solutions in the Non-Linear 
Chiral SU(2) Sigma Model whose original Lagrangian is 



L 1 = ^Jd 3 rTr(d,Ud»U + ) , (1) 

where U is a SU(2) matrix and F n is the pion decay constant. According to 
Derrick theorem || it is necessary to add to expression ([!]) an adhoc stabilizer 
term 
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L 2 = r / d'rTr 

32e 2 



u + d, t u, u + d u u 



where e is a dimensionless parameter. 

The physical properties are then calculated making use of the semi- 
classical expansion of the Quantum Hamiltonian where we perform a rota- 
cional collective coordinate expansion||] U(r,t) = A(t)Uo(r)A + (t), where A 
is a SU(2) matrix and U a static soliton solution. Adopting the hedgehog 
ansatz U = expir.fF(r) where r is the Pauli matrix and F(r) is called the 
chiral angle, the Hamiltonian with the rotational mode can be written as 

H = M+ l -V±V, J = l,2,..., (3) 

where M is the classical energy and / is the inertia moment 0. 

Dube and Marleau [Q indicate therefore a possible way of generalizing 
the Skyrme Model introducing higher-order terms in the derivatives of the 
pion field. In order to simplify the model, they made a particular choice of the 
Hamiltonian parameters that permited to sum the series in an exponential 
funcional form. Using this method they obtained improved numerical results. 

In the present note we propose to introduce another four derivative sta- 
bilizer term given by 

L 2 = J d 3 rc 2 [tr (d^Ud^U^ . (4) 

As this term is the square of the kinetic term defined in (HD, it is natural 
to take it as a pattern for the inclusion of the higher-order derivative terms. 
Thus, the standard form of the Lagrangian terms is 



|rf 3 rc n [Tr(^f/^f/ + )] n , (5) 
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where n=l,2,. ... It must be noted that this form is also invariant under 
chiral transformation (U transforming under U — > AUB^ 1 , A and B are 
SU(2) matrices), and does not destabilize the soliton solution, since we can 
to define a positive static Hamiltonian. 

To simplify the calculations we will adopt the Sugawara form, L M = 
U + d^U, whose the static component can be written as Lj = iT a L®. Using 
the hedgehog ansatz, the kinetic static term (|l|) becomes then, 



Li = — Ci / d 3 rTr 



diUdtlP 



Ci I d 3 r tr [LiLi 



-2ci 



cPrLfLf, (6) 



where L?L° 



2 sin 2 F 



+ F' 2 , and c\ = F%/16. Since the classical energy is 
given by -L, where L is the generalized static Lagrangian, we can write the 
classical total positive energy with the inclusion of the higher order derivative 
terms as 



a Ta\ n 



E = Jd 3 r [ci (2L?L?) + c 2 (2LfLf) + . . . + c n {2L\L 
Jd 3 r 2L?Lf ci [l + a ( 2 L«L») + a (2L"Lf) 2 + . . . + C M2L\L\ 



. n—l 



•(7) 



There are many particular values of the relations — infl7T) that allow the 



energy series to summed. Following Dube and Marleau, if we naively set 



K = — 



\Tl — l ' 



ci ^ri-1)! (2e 2 F2)' 
the series of the classical energy(|7]) converges in an exponential form 



(8) 



E = J d 3 r 2c x Lm exp 



F 1 
16 



Jd 3 r2 



2 sin 2 F 



+ F n 



e 2 F 2 



2 fo dxX 



2sin 2 F 



+ F 
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exp 
exp 



+F r. 



e 2 F 2 
2sm 2 F 



?I2 



(9) 

with i 7 ^ and e being the only input parameters. In expression (|9]) we have 
used the series representation exp y = Y^kLo i 7 ! an d a dimensionless variable 
x defined by 
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x = eF n r. (10) 

It is interesting to point out that with the choice of standard form([5[) and 
with the definition fliTf), all the coupling constants are absorbed in the new 
dimensionless variable x. From (|9|) the variational equation is 



[2x 2 + 2x 2 S + 8x 2 F' 2 + Ax 2 SF' 2 } F" + 8 sin IFF' 2 
+4Ssin2FF' 2 + AxF' - 16Wff ' + AxSF' 

X 

%Ssin 2 FF' 



2sin2F - 2Ssin2F = 



(11) 



where S 



2 sin 2 F 



+ F 



n 



The soliton solution with the baryon number 1 
has the boundary conditions F(x) = it at x = and F(x) = at x — > oo. 
Then, we impose the asymptotic behaviour of F given byfllll) 



lim Fix) 

x— >oo 



B 

~~2 ' 



X 



(12) 



where B is a constant, to obtain using numerical integration the soliton 
solution. In figure 1, we show the numerical behaviour of parameter B, which 
is direclty proportional to the axial vector constant coupling,^ = . 

The inertia moment is calculated similarly to the series of the classi- 
cal energy. Performing the rotacional collective coordinate expansion of the 
classical Lagrangian and picking up only terms linear in Tr(doA + doA) , we 
obtain the expression of the series of the inertia moment given by 



d 3 r- sin 2 F c x [l + 2K 2 (2L^) + . . . + nK n (2L^L^ 
3 



,n-l 



• (13) 



With the substitution of K n by the expression(||), using the dimensionless 
variable x defined in fllOD and taking the limit n oo, we find 
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2tt 1 
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dxx 2 sin 2 F 1 + F' 2 + 



2 sin 



2 FA 



I 



3 F n e 3 



exp F' 2 + 




(14) 




cleon (Mjv = 939Met>) and of the Delta (Ma = 1232Afet>) as input parame- 
ters to determine the pion decay constant F n and the dimensionless Skyrme 
parameter e. The main physical results are shown in table 1, according to[^]. 



Fig 1. Behaviour of the parameter B defmided by B = x 2 F(x) , where F(x) 
is the numerical variational solution of the classical Hamiltonian including 
terms up to order 4(n=2), 6(n=3), 8(n=4) (solid line) and all orders (dashed 
line) in derivatives of pion field. 
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TABLE 1- Physical parameters in the Skyrme Model 





Adkins[|] 


Marleau@] 


This model 


Expt. 


F^(Mev) 


129 


146 


144 


186 


e 


5.45 


8.69 


6.82 




i 

< r 2 >| =0 {fm) 


0.59 


0.60 


0.61 


0.72 


Hp 


1.87 


1.89 


1.90 


2.79 


/'■» 


-1.33 


-1.32 


-1.31 


-1.91 


9 A 


0.61 


0.71 


0.80 


1.23 



Our results indicate an improvement in the physical values of the pion decay 
constant,^, and the axial coupling constant, qa- The others values, i.e., 
the magnetic moments and the isoescalar mean square radius remain basi- 
caly the same obtained by Adkins|| et al. and Marleau[Q. This procedure, 
without doubt, improves the physical results. In order to obtain better re- 
sults for quantities like magnetic moments and the isoescalar mean square 
radius, we should treat the quantization of the classical Hamiltonian (collec- 
tive coordinate quantization in more formal way, using the information 
about constraint that is present in the system 0. This particular study will 
be object of a forthcoming paper 0. 

I would like to thank M.G. do Amaral for critical reading. The work is 
supported by CNPq, Brasilian Research Council. 
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